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Abstract 

Linear and nonlinear Hodge-like systems for 1-forms are studied, with 
an assumption equivalent to complete integrability substituted for the 
requirement of closure under exterior differentiation. The systems are 
placed in a variational context and properties of critical points investi- 
gated. Certain standard choices of energy density are related by Backlund 
transformations which employ basic properties of the Hodge involution. 
These Hodge-Backlund transformations yield invariant forms of classical 
Backlund transformations that arise in diverse contexts. Some extensions 
to higher-degree forms are indicated. MSC2000: 35A15, 58J72, 58A14. 
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1 Introduction 

The study of vectors which are both divergence-free and curl-free can be traced 
back at least to Hclmholtz's analysis of vortices and gradients [14]. The general- 
ization to differential forms which are both closed and co-closed under exterior 
differentiation is the content of the Hodge equations; see, e.g., [35], Ch. 7. The 
divergence-free condition is frequently relaxed in variational contexts, but gen- 
eralizations of the curl-free condition remain rather rare. 

Our goal is to study both linear and nonlinear variants of the Hodge equa- 
tions for differential forms which are neither co-closed nor closed, but which 
satisfy milder conditions having physical and geometric significance. 
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1.1 Organization of the paper 



Sections 1-4 are mainly expository. We introduce the topic in Sec. 1.2 with an 
example from fluid dynamics. Section 2 presents the equations in an invariant 
context. The linear case is studied in Sec. 3, largely as motivation for the con- 
siderably more complex nonlinear case. Two geometric analogies are discussed 
in Sec. 4. Technical results on the properties of solutions are presented in Sec. 5. 
Section 6 shows that ideas introduced, in very different contexts, by Yang [33] 
and by Magnanini and Talenti [19j-|21j. can be given a unified interpretation in 
terms of equations studied in the preceding sections. 

The proofs of Theorems 5, 6, and 8, which are based on rather straight- 
forward applications of nonlinear elliptic theory, are collected in an appendix, 
Sec. 7. We note that these applications are only straightforward once solutions 
have been associated to a uniformly sub-elliptic operator; this is accomplished 
in Sec. 5.2. The methods used in that section to derive uniform estimates are 
elementary and, in particular, do not require a delicate limiting argument that 
has become known as Shiffman regularization (c.f. [33J and the appendix by R. 
J. Sibner to [35]). 

1.2 A motivating example: steady, ideal flow 

In models for the steady, adiabatic and isentropic flow of an ideal fluid, conser- 
vation of mass is represented by the continuity equation 



where v denotes flow velocity and p denotes mass density. The dependence of p 
on |v| 2 is a consequence of compressibility; in the incompressible limit, eq. flj 
says merely that the vector v has zero divergence. If the fluid is irrotational, 
then the velocity is curl-free in the sense that 



Condition @ implies, by the Poincare Lemma, that there exists locally a scalar 
flow potential ip (x) , where x £ M 3 denotes the position of a particle in the flow. 

Perhaps the mildest weakening of the irrotationality condition results from 
replacing ([2]) by the integrability condition 



The replacement of the linear condition ([2]) by the nonlinear condition ([3]) 
as a side condition to eq. flj is likely to result in singular solutions, even in the 
subsonic regime. The usual arguments for reducing p to the conventional form, 
which depend on smoothness (see, e.g., Ch. 1.2 of [4]), would not necessarily 
apply in such cases. This suggests that we consider whether a useful a priori 
bound can be placed on the size of the singular set for solutions of systems 
having the general form {T]), ([3]). In Sec. 5.2 we take the first step toward an 
answer to this difficult question, deriving sufficient conditions under which a 
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V x v = 0. 
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v • V x v = 0. 
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solution remains bounded on an apparent singular set of given codimension. 
Despite the physical motivation (here and in various other examples scattered 
throughout the text), our main interest in this paper is in deriving hypotheses 
which are mathematically natural and apply to large classes of mass densities. 

2 An invariant formulation 

Thus we generalize the mathematical context of Sec. 1.2. Let f2 be an open, 
finite domain of R™, n > 2, satisfying an interior sphere condition. Consider the 
system ([26], Sec. VI; [27], Sec. 4) 

S(p(Q)co) = 0, (4) 



duj = T A lo, (5) 

for scalar-valued 1-forms lo and T, where V is given and lo is unknown; d is the 
flat exterior derivative with formal adjoint 6; Q is a quadratic form in lo given 
by 

Q(uj) = * (w A *w) = {to, to), (6) 

where *: A k — > A"~ fc is the Hodge involution; p is a positive, continuously 
differentiable function of Q (but a possibly singular function of x) . 
Using ([5]), we find that 

lo A duj = -r A lo A lo = 0. (7) 

A 1-form lo that is the pointwise Ricmannian inner product with a vector field 
v is said to be dual to v. In this case the left-hand side of ([7]) is equivalent to 
the left-hand side of ([3]) , and any solution of lo of eq. ([5]) is dual to a solution v 
of eq. ©. 

The system Q, ([5]) is uniformly elliptic provided the differential inequality 

(d/dQ) [Q P 2 (Q)] 

< Ki < j— < n 2 < oo (8) 

P(Q) 

is satisfied for constants Ki,/t2- In the context of fluid dynamics, one typically 
encounters the weaker condition 

0<p 2 (Q) + 2Q P '(Q)p(Q). (9) 

Ideal flow governed by eq. |T]) is subsonic provided © is satisfied. Moreover, 
there is typically a critical value Q cr n such that the right-hand side of © tends 
to zero in the limit as Q tends to Q cr u- In this case, eqs. (01, © with p satisfying 
^ are elliptic, but not uniformly so, and this condition has mathematical as 
well as physical interest. 

If r = 0, then the system ((4]), © degenerates to the nonlinear Hodge equa- 
tions introduced in [36] on the basis of a conjecture in [3]. In that case condition 
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([5]) generates a cohomology class, which is not true in the more general case 
studied here. 

We show in Sec. 2.1 that whenever w is a 1-form, eq. possesses solutions 
of the form 

uj = e^du, (10) 

where u and 77 are 0-forms, and T can be made exact: we can write T = dn. 
When uj is a fc-form, a more general representation applies; that representation 
is discussed in Sec. 2.2. 

In Sec. 5.1 we introduce a variant of eq. (H]) for differential fc-forms (k > 1) 
satisfying eq. with V exact (r = drj), namely 

S [p(Q)uj] = (-l)»( fe +D * (d v A *p(Q)uj) , (11) 

which arises as a variational equation of the nonlinear Hodge energy 

E= \ ( [ Q p(s)dsdfl. (12) 
2 Jo, Jo 



2.1 The Frobenius Theorem 

Let T be fixed, and define 

S= {u e A(fi) = 0A fc (^) : doj = T/\uj}. 
fe=i 

Denote by / = I(S) the ideal generated by S. If T ^ 0, then clearly dl ^ {0}, 
so the fc-forms uj G S do not generate cohomology classes. 

Nevertheless, the ideal / is closed, i.e., dl C /. In fact, a differential form 
a £ / is a linear combination of forms of type uj A (3 with uj £ S, (3 £ A(f2). The 
latter satisfy 

d(ujAf3) = ±uj A (r AP±d0), 

and thus satisfy da G /. This is an important fact, especially for exterior 
systems of 1-forms. 

Following the approach in Sec. 4-2 of [ID] , we define an exterior system {u> a }, 
a = 1, . . . , r of r 1-forms in a space of dimension n = r + s to be completely 
integrable if and only if there exist r independent functions g a , a = 1 , . . . , r such 
that each of the 1-forms uj a vanishes on the r-parameter family of s-dimensional 
hypersurfaces {g a = k a , a = l,...r} generated by letting the constants k a 
range over all r-tuples of real numbers. 

Equivalently, we define {w a }^ to be completely integrable if and only if there 
exists a nonsingular r x r matrix of functions and r independent functions 
g b such that 

r 

uj a = J2^dg b - 

b=l 
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The Frobenius Theorem asserts that an exterior system {w a }^ of 1-forms is 
completely integrable if and only if it generates a closed ideal of A(f2). 

Because a 1-form lo satisfying © generates a closed ideal, by the Frobenius 
Theorem it can always be written in the form (|10[) . (In this case, r = 1.) Thus 

du> = drj A lo, (13) 

which shows that T can be chosen to be exact. Notice the gauge invariance 
having the form T — ► T = T + /(x) lo. 

For this reason, we call the system ([5]) the nonlinear Hodge- Frobenius 
equations for 1-forms. 

Unfortunately, the Frobenius Theorem does not generalize to forms of arbi- 
trary degree k, as the condition dl G I does not imply complete integrability 
if k 7^ 1. However, this does not mean that there is nothing to be said about 
higher-degree forms. Relevant properties of such forms are described in the 
following section. 

2.2 Recursive forms 

An exterior differential form of degree k is said to be recursive with coefficient 
r if it satisfies ©■ Let f2 be star-shaped. It is known that one can define a 
homotopy operator 7i : A k (il) — » A k ~ 1 (fl), which satisfies 

lo = dUoo + Hdw. (14) 

This property can be used, among other things, to show that a closed form on 
a star-shaped domain is exact. We omit the formal definition of this operator 
(for this and further details see Sec. 5-3 of [10]), and only describe its main 
properties: 

(a) 7i is linear 

(b) H 2 = 

(c) HdH = H, dUd = d 

(d) (dUf = dU, (Hdf = Hd. 

Using (b), we observe that Tidui £ KerH. This and (|14|) can be used to define the 
exact part of u) as Lo e = dTiu and the anti-exact part of uj as uj a = cu — uj e = TLdto. 

Using (jT4j) , one can further show that Tt improves regularity. With no loss of 
generality, we prove this for forms having vanishing anti-exact part, i.e., forms 
lo such that 

lo = LO e . (15) 

Note that for any given form lo, no cancelations can occur between co e and u) a , 
therefore Lo e is always as smooth as lo, and furthermore H.Lo a = 0. Thus we can 
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restrict our attention to forms satisfying (fT"5)) . If uj is a one- form, then Ttui is a 
function and (|15p implies that 

= a ■ v«, 

thus improving regularity. In the general case, Tito = (T-Cuj)jdx 1 , where I is a 
multi-index satisfying |7| = k — 1. In order to improve regularity in this case 
one would need to control 6(7iu>), i.e., all derivatives of type d(Hu>) j / 'dx 1 for 
i € I. But S(Hoj) — from the Hodge Decomposition Theorem [25], as Hoj is 
anti-exact (using (6), above); so W is a smoothing operator on fc-forms. 
The important result for us is the following: 

Recursive fc-forms with coefficient T on a star-shaped region can be repre- 
sented as follows [TO] : 

uj = e r > [du + H(9 A du)} , (16) 

where r\ = 9 — HdT, and u — 7i(e _, 'cj). Condition ([5]) implies that 9 
satisfies 

d6 A [du + H(6 A du)] = 0. 



For our purposes we can rewrite (|16j) as 

« = e"5(d«), (17) 

where g is a smooth linear operator, or alternatively, as 

uj = e ri du + e ri h{du) 1 (18) 

where h is also a smooth linear operator, the coefficients of which depend on V; 
the latter variant yields better regularity. In fact, the form h(du) is as smooth 
as u, provided T is smooth. 

A particular case of the above occurs when the coefficient T is exact. In 
that case, the form u> satisfying ([S]) is said to be gradient recursive. For gradient 
recursive fc-forms, (fT6|) assumes the simpler form (p~0|) . in which g is the identity. 



3 The linear case 

Corresponding to the physical example of Sec. 1.2, which illustrates the non- 
linear form of eqs. (J4j) , (J5j) , we can illustrate the linear case by an even simpler 
physical example. 

Condition ([3]) arises when a rigid body rotates in the xy-plane at constant 
angular velocity Cj. Taking the axis of rotation to lie at the origin of coordinates, 
we write the tangential velocity vector in the form 

v = vii + v 2 j, 

where 

Vi = —luu, V2 = Cox. (19) 
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Then 

V x v = 2u)k, (20) 

so © is satisfied; see, e.g., Exercise 4.4 of [24 . In the sequel we take ui = 1 for 
simplicity. 

Equations ([IT))) imply that 

V • v = 0, 

so we express the 1-form u> dual to v as a solution of the linear Hodge-Frobenius 
equations 

Slo = 0, 

duj=TAuj. (21) 

Applying ([T0|) and ([T3)) , we choose r\ to depend only on the distance r from the 
axis of rotation. Then 

dij A uj = rj' (r) ■ r dxdy. (22) 

In addition, eq. ([20]) implies that 

dr\ A u> — duj = 2 dxdy. (23) 

Equating the right-hand sides of eqs. ([22]) and ([23]) . we conclude that ry(r) = 
2 logr and to = r 2 du for u(x, y) = arctan (y/x) . Then \du\ — r _1 , so the singular 
structure of u in the ly-plane is analogous to that of the fundamental solution 
of Laplace's equation in K 3 . In particular, u is singular at the origin of the disc. 

In this example, the Hodge-Frobenius equations themselves are only defined 
on the punctured disc, as 

T Auj = (T\uj2 — T2UJ1) dxdy = duj = Idxdy. 

Using ([ 19[) . we can write this condition as an equation for the inner product 

(&)•(:)-* 

which cannot be satisfied at the origin. 

Thus singular solutions arise naturally in both the linear and nonlinear 
Hodge-Frobenius equations. 

Because eq. (|16[) requires that the domain be star-shaped, the conclusion 
that uj is representable as a product / du, where / is nonvanishing, does not 
follow and is in fact violated in our example, in which / = r 2 . 

We have presented a particularly simple model, in which the role of condi- 
tion (O is especially transparent. For more sophisticated completely integrable 
models of rigid-body rotation see, e.g., [5] and references therein. 

In the linear case we can accomplish easily what we cannot accomplish at all 
in the nonlinear case: an integrability condition sufficient to imply the smooth- 
ness of weak solutions. 
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Proposition 1. Let u> be a weak solution of the linear Hodge-Frobenius 
equations \21\) on fi. If \T\ is bounded and to € L P (Q) for p > n, then u> is 
continuous. 

Proof. The Friedrichs mollification Uh of w is a classical solution of (f2Tj) ; see, 
e.g., Sec. 7.2 of [T3]. Thus 

\<L> h \* + \5u h \* + \w h \* < (\r\p + 

Integrate and apply the i p Gaffney-Garding inequality (Lemma 4.7 of [16]) to 
obtain 

\\Vco h \\ p <C(T)\\iJ h \\ p . 

Because u> £ L p (n), we can allow the mollification parameter h to tend to zero. 
The proof is completed by the Sobolev Embedding Theorem. 

4 Geometric analogies 

4.1 Hypersurface-orthogonal vector fields 

A unification of the two superficially different physical examples of Sees. 1 and 
3 can be found in their underlying geometry — in particular, in the relation of 
each vector field to hypersurfaces created by the level sets of an associated scalar 
function u. We can think of u informally as the potential for the conservative field 
that would result from taking T to be zero in eq. ([5]). Whereas a conservative 
vector field is actually equal to Vm, the vector fields in Sees. 1 and 3 merely 
point in the same direction as Vw. 

A nonvanishing vector field v is said to be hypersurface-orthogonal whenever 
there exists a foliation of hypersurfaces orthogonal to v. The foliated hypersur- 
faces can be represented as level sets of a scalar function it. That is, one can 
write 

v = A(x)Vu 

for a nonvanishing function A. Conversely, a vector field which can be written 
in this way is clearly hypersurface-orthogonal. We conclude that a vector field 
v is hypersurface-orthogonal if and only if the 1-form uj dual to v satisfies 

u> = Xdu , 

with nonvanishing A (i.e., if and only if to is completely integrable). Other 
equivalent conditions now follow from the Frobenius Theorem. 

Hypersurface-orthogonal vector fields arise naturally in general relativity, 
particularly in connection with black- hole mechanics. Introducing a tensor field 

where the semi-colon denotes covariant differentiation with respect to the space- 
time metric connection, the condition that v be hypersurface-orthogonal implies 
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that the antisymmetric part of B, called the rotation tensor, vanishes ([29 , Sees. 
2.32, 2.33). For this reason, vector fields satisfying ([5]) are called rotation-free in 
general relativity, which is somewhat confusing in the context of the examples 
in Sees. 1 and 3. 

4.2 Twisted Born-Infeld equations 

A condition broadly analogous to (O arises if w is a Lie- algebra- valued 2-form 
satisfying the second Bianchi identity. In that case, replacing u> by F A , where 
A is a Lie-algebra- valued 1-form, we have 

dF A =-[A,F A ], (24) 

where [, ] denotes the Lie bracket, an equation which resembles ([5]). 

Precisely, let X be a vector bundle over a smooth, finite, oriented, n- 
dimensional Riemannian manifold M. Suppose that X has compact structure 
group G C SO{m). Let A e F (M,adX <g> T*M) be a connection 1-form on X 
having curvature 2-form 

F A = dA + ~ [A, A] = dA + A A A, 

where [ , ] is the bracket of the Lie algebra 9, the fiber of the adjoint bundle adX. 
Sections of the automorphism bundle Aut X, called gauge transformations, act 
tensorially on F A but affinely on A; see, e.g., 22J. 

Consider energy functional having the form (fT2"| . where Q = |-Fa| 2 = 
(F A ,F A ) is an inner product on the fibers of the bundle adX ® A 2 (T *M) . 
The inner product on ad X is induced by the normalized trace inner product on 
SO(m) and that on A 2 (T * M) , by the exterior product * (F A A *F A ) . 

A nonabelian variational problem analogous to eqs. ([5]), is described 
briefly in Sec. 5.1 of 28J. One is led to consider smooth variations taken in the 
infinitesimal deformation space of the connection and having the explicit form 

var{E)= f p(Q)var(Q)dM = f p(Q)^- \F A+t ^\ 2 dM 
Jm Jm at\t=o 

= f P(Q)i \F A + tD A yj + t 2 ^/\^\ 2 dM, (25) 
J M dt\t=o' 

where D A = d + [A, ] is the exterior covariant derivative in the bundle. The 
Euler-Lagrange equations are 

S (p(Q)F A ) = -*[A, *p(Q)F A ] . (26) 

In addition, we have the Bianchi identity (I24|) . 
Writing eq. {5} in components 

du a =T a b l\uj a , 
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we observe that if — r is interpreted as a connection 1-form, then ([5]) can be 
interpreted as the vanishing of an exterior covariant derivative, which is the 
content of eq. f!?4"]) . Moreover, the well known algebraic requirement that T 
must satisfy 

(dTg -T° A Y%) A w b = 

(c.f. eq. (4-2.3) of [10]) is the zero-curvature condition [Fr,w] = 0. 

If G is abelian, then the Lie bracket vanishes and eqs. (|26|) reduce to the 
system 

5{p[Q(F A )]F A } = 6{p[Q(dA)]dA}=0, 

a nonlinear Hodge equation analogous to taking r = in eqs. Q, (jSJ). 

Equations (|24[) reduce in the abelian case to the equations for the equality 
of mixed partial derivatives, 

d 2 A = 0. 

If p = 1, then eqs. ([26]) are the Yang-Mills equations, describing quantum 
fields in the classical limit. These resemble a version of eq. ([TT]) for 2- forms with 
p = 1, with the Bianchi identity ([24]) playing the role of eq. ([5]). 

Whereas the Yang-Mills equations do not have the nonlinear structure of 
(fTTj) for non-constant p, those of the Born-Infeld model for electromagnetism 
are equivalent to (TTT]) for differential forms of degree 2 with 

P (Q) = {l + \F A f)- 1/2 . (27) 

This model was introduced in [5] in order to produce a model of electromag- 
netism that does not diverge when the source is a point charge. Geometric 
aspects of the model are investigated in [11] and its analytic aspects, in [43]. A 
mathematical generalization of the Born-Infeld model to nonabelian variational 
equations was proposed in [5S] and further studied in [35]; see also [TS] for a 
related problem. The equations assume the form (f2"4"]l . ([2"6"]) for an appropriate 
choice of p(Q). 

Geometrically, the nonabelian model puts a twist in the principal bundle 
corresponding to the configuration space of solutions. Thus eqs. ([24]), ([26]) are 
called twisted nonlinear Hodge equations |28| . A different approach to general- 
izing nonlinear Hodge theory to bundle-valued connections, which is based on 
the formulation of a natural class of boundary-value problems, is introduced in 
[23] ; but interior estimates would be required in order to extend the theory of 
[23] to twisted forms of the equations considered in this paper. The derivation 
of such estimates is a goal of Sec. 5.2. 

5 Analysis 

We do not expect rotational fields of any kind to be very smooth. In particular, 
assumption ([5]) may produce caustics; see, for example, the discussion in Sec. 
2.4 of [29]. However, it is reasonable to seek conditions under which the field 



10 



remains bounded at a singularity or under which the field equations remain 
uniformly elliptic. 

Because these conditions will be derived for a large class of mass densities, 
the strength of the estimates obtained will depend on the integrability with 
respect to x of a given choice of mass density p (Q(x)) . 

In the sequel we denote by C generic positive constants, the value of which 
may change from line to line. We follow an analogous convention for continually 
updated small positive constants e. Repeated indices are to be summed from 1 
to n. 



5.1 Variational structure 

The energy E of the field to on f2 is defined by eq. (|12p . where Q — Q(u>) is 
defined as in © for to E A 1 given by (fTU)) . Then Vi/> S Cfi°(£l), the variations of 
E are computed as 

d 1 f d 

™r(E) = —E(u + ty) |t=0 - - J P(Q)j t Q (« + ^)| t =o dto 

= \ |^(g)|[e 2 "|d( U + ^)| 2 ] |t=0 ^ = 
\ I P(Qy v ^-A\du\ 2 + 2(du,tdi/j)+t 2 \diP\ 2 ) dn = 



2Jn yJ dt Vl 1 x ' r/ 1 rl y l*=0 
p{Q)e 271 {du,dip)dn= ( (p(Q)e 2v du,diP)dn 



d{p{Q)e 2,1 du^)dVL+ / (<S [p(Q)e"w] , ^)dVL = / ((5 [p(Q)e n Lj] ,^)dSl, 

as ip has compact support in Q. At a critical point, var(E) = 0, or 

S [p(Q)e r 'to} = 0, (28) 

so the variational formulation yields a "weighted" form of the continuity equa- 
tion ([?]). The presence of this weight adds an inhomogeneous term to the un- 
weighted variant. To see this, we write the local equation 

-d t [p(Q)eV] = - (dm) e"p(Q)w< - ^d, \p{Q)ui\ = 0. 

This local form corresponds to the invariant representation 

e r >5[p(Q)to} = e*>p(Q)(dr l ,io). 

Thus we obtain eq. for one- forms to (i.e., for k — 1) or equivalcntly, since 
(5ce = — * d * a, Va G A 1 , the equation 

d * [p(Q)to] = -drj A * [p(Q)to] . (29) 
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5.1.1 Variational equations for A;- forms 

Recall that an exterior differential form u> of degree k which is recursive with 
coefficient T can be written in the form (fT5|) , where u is a (k — l)-form (which 
depends on lo), and where the function 77 and the linear operator h depend only 
on T (Sec. 2.2). 

We now compute the variation of the energy functional (fT2|) among all forms 
u + ta satisfying 

d (lj + ta) = T A (oj + ta) . 

Such forms satisfy (fT5|) with fixed 77 and h. Therefore the variations of E are 
computed as 

1 \ f d 

var{E) = —E(uj + ta) ]t=Q = - p(Q)-fiQ (w + ta) [t=0 dfi 

(p(Q)uj,a)dn, 



n 



where a = e 1 ' [dv + h(dv)) for a (fc — l)-form v which depends on a. 
Thus 

var(E) = - [ (e n p(Q)u;,dv + h(dv))dn = 
2 Jn 

= [ (e ri p(Q) U j,Gdv)dn= [ (G T e r 'p(Q)uj,dv)dn, 
Jn Jn 

where G = ffij(r) is an n x n matrix and G T is its transpose. 

As the forms a (and thus v) are assumed to have compact support in J7, 
setting var(E) — is equivalent to imposing the condition 

6 [G T e n p{Q)u] = 0. 

Notice that if lj is gradient recursive, then G is the identity matrix and we 
recover (|2"gj) : or equivalently, since Va G A fe 

6a = {-l) nk+n+1 *d*a, (30) 

the equation 

d * [p{Q)uj] = -dr) A * [p(Q)to] , (31) 

i.e., eq. (jlip for a gradient-recursive fc-form u with coefficient £^77. 

The variations employed in this section, applied directly to a 2-form, are 
necessarily different from the variations of (|25p. which are applied instead to 
the Lie-algebra valued connection 1-form A. 

In the remainder of the paper we focus mainly on 1-forms. However, many 
of the results will extend easily to gradient-recursive fc-forms. Moreover, the 
structure of expressions fTT)) and (fTB)) suggest that, under appropriate technical 
hypotheses on the linear operators g or h, many of the results will extend also 
to recursive fc-forms with coefficient T (not necessarily exact). 
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5.2 When are solutions bounded at a singularity? 

Following [32]) we find it convenient to introduce a function H(Q) which is 
defined so that 

H'(Q) = \p{Q) + Qp'{Q). (32) 

Then ellipticity is equivalent to the condition that H has positive derivative 
with respect to Q. 

In Theorem 7 and Corollary 8 of [26], and in Theorem 6 and Corollary 7 
of [57], L p conditions are derived which imply the boundedness of solutions to 
eqs. dU, on domains that include singular sets of given co-dimension. We 
call such theorems partial removable singularities theorems, as they imply that, 
although solutions may have jump discontinuities at the singularity, they cannot 
blow up there. Those results require the mass density p to satisfy the inequality 

C{K + Q) q <H'(Q)<C- 1 {Q + K) q (33) 

for constants q > and K > 0. (This hypothesis is imposed in [27J, following 
Sec. 1 of [42]; a somewhat stronger hypothesis is imposed in [26].) 

There is obvious interest in deriving estimates for densities which may not 
satisfy ([33]) . The focus of this section is to obtain and exploit such estimates for 
a broad class of densities, using the variational form of the equations. However, 
we retain the condition that p is positive, which is natural for applications. 

We impose an additional condition that arises from technical considerations. 
If p'{Q) > 0, we require that 

H'(Q) < C P (Q). (34) 

(This inequality is satisfied automatically if p'(Q) < 0.) If p'(Q) < we require 
instead that 

H'(Q) > Cp(Q). (35) 

(This inequality is satisfied automatically if p'(Q) > 0.) Note that (j35|) implies 
([9]) under our assumption p{Q) > 0. In the case p'(Q) < 0, eq. ([TT]) is uniformly 
elliptic (i.e. ([8]) is satisfied) whenever (|35p is satisfied and p is noncavitating: 
bounded below away from zero. In the case p'{Q) > 0, condition ([8]) is satisfied 
whenever p(Q) is bounded above. 

Densities which satisfy ([33]) satisfy the hypotheses ([34]) . ([35]) . However, there 
are many densities which satisfy ([34]) . ([35]) but do not satisfy ([33]) . Among the 
latter are densities for which the value of the exponent q on the left-hand side of 
inequality (|33| differs from its value on the right-hand side, and certain densities 
for which the value of q in ([3"3"|) is negative. 

As a simple illustration, consider the class of densities 

p(Q) = (K + Q) 9 , -1/2 <q<0,K>0. (36) 

Such densities do not satisfy condition ([33]) and cavitate as Q tends to in- 
finity. They arise, for example, in connection with models of pseudo-plastic 
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non-Newtonian fluids [2]. In this section we will obtain sufficient conditions 
for a bound on Qp(Q) which is valid as Q tends to infinity. Under additional 
conditions, this bound extends to possibly singular solutions of the system ([5]), 
(fTTj) with density given by (1551) . Note that for such densities, a bound on the 
product Qp{Q) implies an asymptotic bound on the norm Q of the solution 
itself. 

Because conditions (|34l) and (j35|) are used in a crucial way to establish both 
the subelliptic estimates and the ellipticity of the second-order operator, they 
appear to provide a mathematically natural generalization of condition (|33[) . 

Lemma 2. Let condition \35\) be satisfied. Then H can be chosen so that 

Qp{Q) < CH{Q). (37) 

Proof. If p'(Q) is non-negative, then ([57]) is always satisfied with G = 2 if 
we choose H(0) > 0. (Note that assumption (j3"5|) is not needed in this case.) In 
order to see this, let 

$(Q) = 2H(Q)-Qp(Q). 
Then $(0) = H(0) > 0, and 

$'(Q) = 2H'{Q) - p(Q) - Qp'{Q) = Qp'(Q) > 0. 

Thus $(Q) remains nonnegative on the entire range of Q. 
If p'(Q) < 0, we assume (|35p. Then in particular, 

2ff'(Q) > ep(Q), (38) 

where we take e to be so small that it lies in the interval (0, 1) . Inequality (|38[) 
can be written in the form 

(1 - e) p(Q) + 2Qp'(Q) > 0. (39) 

Define a constant c by the formula 

1 + e 
C= —- 

In terms of c, (|39p can be written in the form 

(c - 1) p{Q) + 2 (c - 1/2) Q//(Q) > 0. 
We can convert this expression into the differential inequality 

2cH'(s) > 4- M*)] > * e [0,Q]. 
as 

Integrate the inequality over s, using H(0) — 0. We obtain l[3"7]). 
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Lemma 3. Let the 1- forms T and u> smoothly satisfy JSJ) and (!!])■ Let 
p > satisfy conditions \3Jj\ ) and H35\) . Then 

AH + (-1) 3 "V • {* [u A * { P '(Q)dQ A u)]} + C (\T\ 2 + |Vr|) H>0. (40) 

Proof. We have (c.f Sec. 1 of 02]) 

(w, A WQH) = d % {uo, d % (p(Q)cu)) - (d.cu, di (p(Q)lo)) 

= AH(Q) - [p(Q)(diUi,diUi) + p'iQWiuMdiQ] > (41) 

where 

AH(Q) = di (diH(Q)) = 8 t (H'(Q)diQ) . 
Writing diQ = 2(diUj, u), we rewrite (l4"Tj) in the form 

(w, A [p(Q)cu]) = AH(Q) ~ p(Q)\Vlu\ 2 - 2Qp'(Q) \d\u\\ 2 . (42) 

Applying eq. (j4"2"|) to the operator identity A = — (Sd + dS) and using (fTTj) , we 
write 

= (ui, A [p{Q)w]) + (u, Sd (p(Q)oj)) + (u, dS (p(Q)ui)) 

= AH(Q) - 1 + {w 1 5 {dp A u)) + n + r 2 , , (43) 

where 

= P {Q)\Vlu\ 2 + 2Qp'{Q)\d\u\\\ 
n = (u,5(p(Q)du)), 

and 

r 2 - («,<i|/,(Q)<I»]). 

Define 

L U {H) = AH + (uj, S (dp{Q) A to)). 
Then f|43[) can be written in the compact form 

LuH + ti + r 2 = 7. (44) 

We have 

n = (w, p(Q)6 (T A w) - (dp(Q),du)) 
< \u>\[p(Q)\5(rAcj)\ + \dp(Q)\\dw\] < 
M b(Q)|VT|M +p(Q)|r||Vo;| + |^(Q)dQ||dw|] . 
Applying (|5|) to the last term on the right-hand side and using 

|dQ|=2M|dM|, 
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we obtain 



ri < Qp(Q)\VT\+ p(Q)\T\\u;\\Vu\ + 2Q|^(Q)||r||o;||d|w|| < 



Qp(Q)\VT\+p(Q) 



Q\£l 

2s 



+2Q\p'{Q)\ ( ^-\d\u;\\ 2 



iVwl 



2e 



(45) 



Estimating T2 yields the same terms: 



T2 < Qp(Q)|Vr| + |«||r| [p(Q)|Vo;| + MVp(Q)|] = 
Qp(Q)|Vr| + |w||r|p(Q)|V«| + 2Q|r||//(Q)|H|d|w||, 

which is bounded by the right-hand side of (|45|) . That is, 

ri + r 2 < 2Qp(Q)|Vr| + 9. [ p (Q) + 2Q\p'(Q)\] |lf + 

£ 

ep(Q)\Vv\ 2 + 2Qe\p'(Q)\\dH\ 2 . 



(46) 



(47) 



Applying inequality (|47|) to eq. (|44|) , we obtain, in the case p'(Q) > 0, the 
estimate 



L U (H) + 2 



Q P (0)|vr| 



QH'(Q). 



> 



(1 - e) [p(Q)|V W | 2 + 2Qp'{Q)\d\u\\ 2 } > (1 - e) i/'(Q)|d|o;|| 2 , 



(48) 



the inequality on the right following from Kato's inequality. We apply Lemma 
2 and (|34ll to terms on the extreme left-hand side of inequality (f48|) : 



Qp(Q)|vr| + ^M|rf 



< 



c[tf(Q)|vr| + Q P (Q)|r| 2 ] < 
c(|vr| + |r| 2 )i?(Q). 

We now have, for the case p'(Q) > 0, the estimate 

L U (H) + c (|vr| + |r| 2 ) H > (1 - e) i?'(g)|d|w| 2 > o. 

In the case yo'(Q) < 0, we also apply (|4T|) to but here we obtain 



(49) 



g P (Q)|vr| + ^M|r| 2 



> 



(1 - e) [p{Q)\Vuo\ 2 + 2Qp'(Q)|d|w|| 2 ] + 4 £ Qp'(Q)|d|a;|| 2 
In the case p'(Q) < 0, condition ([55)) yields 



(50) 
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thus 

Q\p'(Q)\ < CH'(Q). (51) 

(Note that (|5"Tj) is automatic for p'(Q) > 0.) Applying Kato's inequality and 
(ED) to (50} yields 



LUH) 



Qp( Q )\vr\ + 9I£M lTl 2 



> 



(1 -e-AeC) H'(Q)\d\oj\\ 2 > 0. 

Applying and ([57]) to the left-hand side of this inequality we obtain (I49[) 
(for updated e) for the case p'{Q) < as well. 

It remains only show that the operator L w (i?) can be put into divergence 
form, at the cost of absorbing another lower-order term. The second term of 
L UJ {H) can be written in the form 

(oj, 5 {dp(Q) A uj)) = * [uj A *S (dp(Q) A lj)] = 

*d [uj A * (dp(Q) Au>)]-* [duj A * (dp(Q) A uj)} , 

where 

— * [duj A * A w)] = - * [r A uj A * (dp(Q) A uj)] > 

-2Q|r||p'(Q)|M|d|u;||, 

which is estimated in the same way as the last term in the sum on the extreme 
right-hand side of (|46|) . 

Taking into account (|30[) we obtain, for any /c-form a, 



*da = (-i)Hn-k) * d (^ a = 
(-l) fe ("- fc ) (*d*) * a = (_ 1 )2fc"+™+i-fe 2 j * a (52) 

Taking fc = 1 and 

a = *[uj A * (p'(Q)dQ A uj)] 

in eq. I|52p , we can express the operator S in that equation as a divergence. This 
allows us to write P0|) (again updating C and e) in the form 

AH + (-l) 3n V ■ {* [u A * (p'{Q)dQ A w)]} + C (|rf + |vr|) H > 

(1 -e) ff'(Q)MM| 2 > 0. 
This completes the proof of Lemma 3. 

Lemma 4. Under the hypotheses of Lemma 3, the operator 
C U (H) = AH + (-1) 3 "V • {* [uj A * (p'(Q)dQ A uj)}} 
is a uniformly elliptic operator on H. 
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Proof. We argue as in Sec. 1 of [12] and Sec. 4 of [57], but without using 
hypothesis ([33]) . 

Define a map /3 U : A — > A fc+1 by the explicit formula 

Aii : A* ~~ * A w, (53) 

for /i£ A" and w S A fc . Then we can write the variational form of the Hodge- 
Frobenius equations {5J, ([TU]) . ([TT]) in the alternate form 

dxTi = -p ai (rj), a = 1,2, (54) 

where o~\ — *p(Q)ui and 1T2 = — The "irrotational" case <iw = can be 
recovered as the special case of ([5"4"]) in which 

A«(»?) =d(r)uj) , 

in which case eq. (|53[) implies that 

dr] A u> = drj A ui -\- r/dio 

(c.f. Sec. 4 of [27] and Sec. 1 of 02]). 
Moreover, writing 

Aw (5) = dg Auj 

for some 0-form g, we compute for arbitrary compactly supported p, € A fe+1 , 

(//, dg Auj) = * (dg A (**) (w A = (dg, * {lo A *p)) = (g, S * [to A *//)). 
So the map ^* : A k+1 — ► A defined by the explicit formula 

Aw(M) = <5*("A*At) 

is the formal adjoint of (3 U . 

In terms of the maps j3 u and /3*, we can write 

(-1)" +1 V ■ {*{u A * (p'(Q)dQ A 00)}} = f3*Mp] 

= AtM#)] 

for piu satisfying 

^) = ^ )d HA.. 
Using (3* , we write the inequality of Lemma 3 in the form 
£uj(H) + lower-order terms > 0, 

where 

C u (H) = AH-fc\ Jiu (H)]. 

Writing 

C u (H)=d k (a jk d j )H, 
18 



we find that if p'(Q) < 0, then (|51[) implies that 



y< Q\p\Q)\ , Mr 

1 < a kj H < ajy + G. 



If p'(Q) > 0, then we write 

Qp'(Q) 



v • 



l - 



H'(Q) 



H'(Q) 



vh y = v 



2ff'(Q) 



Vff ; 



condition (|34[) implies that there is a positive constant c such that 

c g(g) 

2 " 2JT'(Q) " ' 

c is the reciprocal of the constant C in (fM|) . This completes the proof of Lemma 
4. 



These three lemmas easily yield: 



Theorem 5. Under the hypotheses of Lemma 3, the product Qp{Q) is locally 
bounded above by the L 2 -norm of H. 



The proof of Theorem 5 is given in Sec. 7.1. 



Of course the integrability of H depends on p. But for any given p in C 1 (Q), 
H can be computed explicitly by integrating (|32p . Even if p cavitates, Theorem 
5 yields asymptotic information about the fastest rate at which Q can blow up. 
Nevertheless, Theorem 5 is ultimately not very useful, due to the hypothesis 
that the solutions are smooth. It would become more useful if it could be ap- 
plied to singular solutions. We will find that the partial removable singularities 
theorems proven in [27] for singular sets of prescribed codimension extend to 
our conditions on p under slightly different hypotheses. 

Initially, we treat the special case of an isolated point singularity, for which 
the proof is somewhat simpler than the proof for higher-order singularities and 
the range of applicable dimensions somewhat larger. 

Theorem 6. Let the hypotheses of Lemma 3 be satisfied on f2\{p}, where p 
is a point ofR n and n > 2. If H £ L 2 "/( n ~ 2 )(f2) and if the function 

/ = |vr| + |r| 2 (55) 

is sufficiently small in L"/ 2 (f2), then H is an H 1,2 -weak solution in a neighbor- 
hood of the singularity. 



The proof of Theorem 6 is given in Sec. 7.2. 
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Corollary 7. Let the hypotheses of Theorem 6 be satisfied and, in addition, 
let the function f given by \55\) satisfy the growth condition 

f \f\ n/2 dn < Cr K (56) 

is r (x )nf2 

for some n > 0, where B r (xo) is an n-disc of radius r, centered at Xo. Then the 
conclusion of Theorem 5 remains valid. 

Proof. Apply Theorem 5.3.1 of |25j to the conclusion of Theorem 6, following 
the proof of Theorem 5. 

Note that the singularity in Corollary 7 is in the solution, rather than in the 
underlying metric (c.f [39| . in which metric point singularities are considered 
in the case T = 0) . Corollary 7 extends to higher-order singularities in spaces of 
sufficiently high dimension; but the proof requires more delicate test functions. 

Theorem 8. Let the pair u), T smoothly satisfy eqs. hi 1\). with p 

satisfying the hypotheses of Lemma 3, on fi\E. Here £ is a compact singu- 
lar set, of dimension < m < n — 4, completely contained in a sufficiently 
small n-disc D which is itself completely contained in the interior of tt. Let 
H(Q) lie in L 2 ^ (D) n L 272 (D), where (5 = (n - m - e) / (n - m - 2 - e) for 
1/2 < 7! < 7 2 . If the function f given by [55]) satisfies the growth condition 
i56\). then the conclusion of Theorem 5 remains valid. 

The proof of Theorem 8 is given in Sec. 7.3. 

The conditions imposed on p in this section also lead to extensions of known 
results for the conventional case T = 0. In particular, we consider equations 
which, expressed in components, have the weak form 

/ [p(Q)u Xk } x ^ Xk dn = 0, (57) 
Jn 

where ip £ C^(fl) is arbitrary and Q = \du\ 2 . Such equations have been in- 
tensively studied in cases for which p [\du\ 2 J grows as a power of \du\; c.f. eq. 
(3.10) of [5]. Equation (f5"T)) can be interpreted as a weak derivative with respect 
to Xi of the system gj, © with r = 0. 
Define 

fi(Q) = Q P 2 (Q). 

The system (T?]), ([5]) is elliptic precisely when Sj'(Q) > 0. The following result 
extends Theorem 1 of [37] , which requires the derivative of p with respect to Q to 
be nonpositive. If p is specified to be nonincreasing in Q, then the integrability 
of $)(Q) can be shown to follow from finite energy. In the general case, we 
impose this integrability as an independent hypothesis on the weak solution. 
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Theorem 9. Let the scalar function u (x) satisfy equation |57[ ) with p 
bounded, positive, and noncavitating, and with Sj(Q) integrable. Assume condi- 
tions (OJ) and \3$ - Then for every n-disc Dr of radius R completely contained 
in f2 there is a positive number 6 > such that 



sup Si(Q) < CR- n { «(Q) * 1, 



where C depends on p and 8 but not on Q or R. 

Proof. Define Q* to be the (possibly infinite) value of Q for which 

sup p(Q) = p (<2») . 
Qen 

Then 

where C is the constant of (|34p. 

As in |27j . we initially estimate smooth solutions, and subsequently extend 
the result to weak solutions, recovering the derivatives as limits of difference 
quotients. 

We choose test functions having the form 

^(x) = u x% h al \\ 

where a > 0, £ (x) S Cg° (D R ) , and 

fi = ft(Q)+e (59) 
for a small positive parameter e (c./. [8], Sec. 3). Then 

[p{Q) u x k ] Xt vi k = 

2H'(Q) (u XtXk f f) a/2 ( 2 + ap(Q){H'(Q)} 2 \VQ\ 2 h {a - 2)/2 ( 2 
+2H'(Q)Q Xk h a/2 (Q k =n+i2+ i 3 . 

Here 

k > 2H'{Q)u XiXk u XtXk Qp 2 {Q)h {a - 2)l2 C, 2 = 
\p 2 {Q)H'{Q)\VQ\ 2 h {a - 2)/2 C, 2 > C\X7 



where the last inequality follows from (|58)l and the constant C depends on a. 
P^ 1 (<?*) , the constant of (|58|) . and the lower bound of p{Q). 
Similarly, 

p{Q) 



«2 > 



p(Q*) 



12 > c|v (h {a+2)/i ) | 2 c 2 
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and 

/:; ■2II'U))Q, -;:X;C, • 

The latter quantity can be estimated by applying, as in the proof of Lemma 3, 
the elementary algebraic inequality (Young's inequality) 

ab > - fa + l6 2 ) , 

now taking 

a = 2H\Q)h {a - 2)/ \Q X] , b = h [a+2)/ \ X] , 

and e = 5p 2 (Q) for small 5 > 0. 

Putting these estimates together, we obtain 

f | V ^(«+2)M |2 ( 2 d0 < C / f,(«+2)/2|VC| 2 ^. 

Jn ^ ' Jn 

The proof for the smooth case is completed by applying the Moser iteration as 
in expressions (9.5.8)-(9.5.12) of [TB] and subsequently letting the parameter e 
in d59"|) tend to zero. The proof is extended to the general case by applying the 
difference-quotient method as in eq. (13) of [57]; see also Lemma 2 of [12"] . 

6 Hodge-Backlund transformations 

Different choices of p may sometimes be related by a special kind of Backlund 
transformation which is based on properties of the Hodge involution. We call 
these transformations Hodge-Backlund. Although this term does not seem to 
have been used up to now, such generalized Backlund transformations have a 
long history in diverse helds of mathematical physics. Our aim in this section is 
to unify these various transformations, place them in an invariant context, and 
extend them to the completely integrable case. 

Historically, the term Backlund transformation has been defined in many 
ways; see, e.g., [31] for the classical theory. In the sequel we will use it in the 
general sense of a function that maps a solution a of a differential equation A 
into a solution b of a differential equation B and vice-versa, where B may equal 
A but b will not equal a. 

6.1 Transformation of the Chaplygin mass density 

The mass density for the adiabatic and isentropic subsonic flow of an ideal fluid 
has the form 

p{Q) = f i - —Q) . (6°) 

for Q £ [0,2/ (7 + 1)) , where 7 is the adiabatic constant: the ratio of specific 
heats for the gas. The adiabatic constant for air is 1.4. Choosing 7 to be 2 
we obtain, by an independent physical argument originally introduced for one 
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space dimension in [30], the mass density for shallow hydrodynamic flow in the 
tranquil regime (c.f. eq. (10.12.5) of [40]). If we choose 7 to be —1 (a physically 
impossible choice) , we obtain the density of the minimal surface equation |17] , 

P{Q) = vTTt? (61) 

Flow governed by this density is called Chaplygin flow. Despite the fact that 
the numbers —1 and 1.4 are not particularly close, this choice of mass density 
has many attractive properties as an approximation for (|60[) ; see, e.g., [9] and 
Ch. 5 of |4J. These properties are, in general, retained in the case of completely 
integrable flow described in Sec. 1.2. 

If T = 0, eqs. (gj and (0 with p(Q) given by (JM]) describe, for k = 1, 
nonparamctric minimal surfaces embedded in Euclidean space. If k = 2 they 
describe electromagnetic fields in the Born-Infeld model, as in (|2"T|) . 

More generally, we have the following result, which extends an argument 
introduced for the case T = by Yang in [43]; see also [1] and Theorem 2.1 of 



Theorem 10. Let the 1-form lj satisfy eqs. f5]) and ill]) , with p satisfying 
H61\). Then there exists an (n — l)-form £ with |£| < 1, satisfying equations 
analogous to (E|) and \11\) . but with V = drj replaced byT = dr) = —drj and p(Q) 
replaced by 

p(\^\ 2 ) = ^^=. (62) 

V 1 - l£l 2 



Proof. Equation (fTTj) can be interpreted as the assertion that the (n — 1)- 
form 



satisfies 



UJ 



(63) 



d£ = df) A £, (64) 

that is, equation ([5]) with r\ replaced by r) = —r\. As a consequence, we conclude 
that this (n — l)-form is also gradient recursive and, on domains with trivial de 
Rham cohomology, there exists an (n — 2)-form a such that £ = e^da (c.f. Sec. 
2.2, final paragraph). 

Because the Hodge involution is an isometry, 

I, ,|2 



1 

or, equivalently, 
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Note that eq. (j6"5)) implies |£| 2 < 1, as well as 



i i 



= 1. 




(-1) 



n-1 



U). 



Hence, 



d*m\ 2 )0 = (-l)"- 1 ^ = (-lf-^Au; = d V A* (/5(|e| 2 )e) = —dfjA* (p(|£| 2 )£) , 



which is equivalent to equation for the gradient recursive (n — l)-form £ 
with coefficient d?), where p has been replaced by p (see also ([31]) ) . This 
completes the proof of Theorem 10. 

The above argument carries over to any pairing of functions /?(|a;| 2 ), p(|£| 2 ), 
as long as their product is 1. Moreover, it carries over essentially unchanged 
to gradient-recursive k- forms to, as these would automatically yield gradient- 
recursive (n — fc)-forms £. 

The same argument, with small modifications, extends Theorem 10 to gen- 
eral (non-gradient-recursive) £;-forms, in which we write equations ([5]) and pip 
in terms of T and T rather than in terms of 77 and r). 

In addition to the original, "irrotational" version of Theorem 10 introduced 
in [43] , other aspects of the duality of mass densities for nonlinear Hodge equa- 
tions are presented in [36], and in Sec. 4.2 of [28]. (The irrotational case of the 
above argument is reviewed in Sec. 4.1 of [28] : note the recurring misprint in the 
two paragraphs following eq. (35) of that reference: dio should be du.) Densities 
of the form (|62p arise in the study of maximal space-like hypersurfaces [7] and, 
in a completely different way, harmonic diffeomorphisms [41] . 

6.2 Transformation of the complex eikonal equation {after 
Magnanini and Talenti) 

As an example of the diverse fields in which these transformations arise, and of 
the simplifying and unifying role of the Hodge involution, we describe an exam- 
ple from complex optics. The description of the local Backlund transformations 
follows the analysis of Magnanini and Talenti [H5] , in which these transforma- 
tions were introduced; see also [20] and [21]. We then reproduce the argument 
of [19] in a simpler, invariant context using the Hodge operator. 
The eikonal equation in K 2 can be written in the form 



where v (x, y) is a given real-valued function. If we write the solution tj) (x, y) 
as a complex function having the form 



ifc + ^ 2 + v 2 = 0, 



(66) 



ip{x,y) 



u (x,y) + iv (x,y) 
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for real- valued functions u and v, then (|66[) is equivalent to the first-order system 

(67) 



2 i 2 2 2 i 2 r\ 

U x +U y -V x -V y +V = 0, 



U X V X + UyVy = 0. (68) 

The function v corresponds physically to the refractive index of the medium 
through which the wavefront represented by the function tp propagates. If 

,2 i „,2 



+ Uy > 0, 



(69) 



then v x and v y cannot both vanish, by eq. (|67|) . In that case, eq. (|68|) can be 
expressed, in the language of proportions, by the assertion that cither 



v l + v l = ± 



Ul+U 2 y, 



(70) 

u x = u y = 0. (71) 

Treating these relations as a coupled system of scalar equations, we have, 
under the same hypothesis, 

= ±- 



\j V l+ V l \j U l + 



and, using eq. (|67] 



U-B = i (~U y ) 



(72) 



Similarly, the scalar equation 



±- 



vl+vl 



\J u l + u l 

which also follows from (|70p . implies by analogous operations the equation 



Vy = ±U X 



Ul+U 2 y 



(73) 



We can write the coupled system (|72")l , (|73p as a vector equation of the form 

(74) 



V x 









Ul+U 2 y 



or as a single (exact) equation for 1-forms, 



du = ±W1 



(—Uydx + u x dy) 



(75) 
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This implies the local existence of a solution to the divergence-form equation 

(76) 



K + u y 



0. 



whenever condition (f6"9"| is satisfied. 

The equations (RTF]) , (f6"8")l . and ([70)1 can also be solved for u x and u y , in 
addition to being solvable for v x and u y as in l|72p and ([75)1 . Under the hypothesis 
that either 

(77) 



2,2 2 
+ V y = V 



or 



2 i 2 ^ 2 
^ + v y > V , 



J x ' "y 

one obtains, by completely analogous arguments to those applied to v x and v. 
the vector equation 



(78) 



= =Fa 1 



(79) 



and the equation 



du = =F 4 / 1 



v 2 4- i> 2 v 9 



(—v y dx + v x dy) 



(80) 



for 1-forms. Note that eq. (|80[) is exact if either (|77|) holds, or if (|78|) holds and 
the divergence-form equation 



dx 



'l - 



1? + 1? 



d_ 

dy 



'1- 



t; 2 + u 2 y 



= 



(81) 



is satisfied. 



6.2.1 Hodge-Backlund interpretation 

Equations (|76|) and (f8Tjl define a Backlund transformation u — > t> and its inverse 
v — ► u (c./. |21jV These equations can be written in the form of nonlinear Hodge 
equations 

d*(p(Q)w) = 0, (82) 
= (83) 
(that is, as eqs. (@J , © for T = 0) . Either 

km 2 ) - yi+^J (84) 
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corresponding to eq. (|?6"j) . or 



mn = ^- w (85) 

corresponding to eq. (|8ip . In either case we assume that Q does not vanish, by 
analogy with the equations J69]), ([71]), ([77]), and (|78f . 

Take p as in (|84|) . If the domain is simply connected, then eq. ([55)) implies, 
analogously to Sec. 5.1, that there is a 0-form r such that 

uj = du, 

and a 0-form v such that 



d« = ±* Wl + T-ijw =±* #f+^ n . (86) 



Because the Hodge operator is an isometry — which is illustrated locally by 
(j70[) , we have 

\dv\ 2 = \du\ 2 + v 2 . 



Thus the Hodge-Backlund transformation (|86l) yields an invariant form of (J67J). 
Unlike classical Backlund transformations of the eikonal equation, in this case 
the Cauchy-Riemann equations are not satisfied. Rather, 

U x = TP(Q)Vy 

and 

Uy = ±p(Q)V X , 

which is sufficient for the orthogonality condition 

Now take p as in ([85]) . Arguing as before, we conclude that there is a 0-form 
u such that uj = du, and a 0-form v such that 



V M 

We obtain 

\dv\ 2 = \dii\ 2 - v 2 . 

Letting u — ±iu and v = ±iv, we obtain a mapping taking solutions to eqs. ([5]), 
(fTTj) with r = and p satisfying ([85]) into solutions of that system with Y = 
and p satisfying ([84]) . 

These arguments extend immediately to the Hodge-Frobenius case — for 
example, by replacing (|86| with the expansion 



e^dv = ± * (^e 2 ^\du\ 2 + v 2 ^j 
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and squaring both sides. They also extend in a straightforward way to gradient- 
recursive /c-forms and, with some modifications, to general fc-forms (see the 
remarks following Theorem 10). 

Motivated by these examples, we offer a general definition of the Hodge- 
Bdcklund transformation. It is a map taking a solution a of a nonlinear Hodge- 
Frobenius equation having mass density pa into a solution & of a nonlinear 
Hodge- Frobenius equation having mass density ps and vice- versa, where B may 
equal A but b will not equal a. 

7 Appendix: Methods from elliptic theory 

In this appendix we collect the proofs of Theorem 5, 6, and 8, which follow 
directly from the association of solutions to a uniformly subclliptic operator via 
Lemmas 2-4. 

7.1 Proof of Theorem 5 

We require the following well known extension of de Giorgi-Nash-Moser theory: 

Theorem A.l (Morrey [25], Theorem 5.3.1). Let n > 2. Let U S H^ 2 (D) 
for each n-disc D CC ft, where f (x) > 1 and define an L 2 -function W = U x 
for some A G [1,2). Suppose that 

[ (a a0 dpWd a ( + fW()dri<Q (87) 
Jn 

y( G Cq°(Q) with C( x ) > 0, where the coefficients a and f are measurable with 
f G L n / 2 (fi). Let the matrix a satisfy the ellipticity condition 

C\i\ 2 < a a "(x)£ a & 

for |a(x)| < M at a.e. x G £1 and for all £. Moreover, let the growth condition 
A 56]) be satisfied. IfUE L 2 (fl), then U is bounded on each n-disc D CC and 
satisfies 

|C/(x)| 2 <Ca » [ \U(y)\ 2 dy, 

where x is a point of Dr (xo) C O. 

Slightly modified conditions will extend Morrey 's result to n = 2; see [2"5] . 
Sec. 5.4. 

We apply Theorem A.l, taking U(x) = iJ(Q(x)) + 1 and A = 1. Let a be 
given by the matrix a of Lemma 4 and let / be given by ([53)1 . Then / satisfies 
the local growth condition (|56J) by smoothness. The inequality of Lemma 3 
completes the proof of Theorem 5. 



2N 



7.2 Proof of Theorem 6 

Without loss of generality, take p to be the origin of coordinates. Write the 
operator in the form 

Co; (H) = AH + V • [T(H)] 

for 

T = (-1) 3 " * {u A * [p'(Q)dQ A lu}} . 

Then 

V [T(H)] = d k (a jk djH) 
by arguments analogous to those of Lemma 4. Moreover, Lemma 4 implies that 

\VH\ 2 < (1 + a jk ) djH ■ d k H < C\VH\ 2 . (88) 

If / is given by then 

/ { [ViJ + T(H)] ■ VC - CJHQ} dn<0 (89) 

for all nonnegative test functions ( 6 C^(£l) which vanish in a neighborhood 
of the origin. Define 

C = (vVuf H i 

where rj G C£° (D) for an n-disc D of radius R, containing the origin, and itself 
completely contained in the interior of fi. Recalling that we have taken H (0) = 
and that ^ implies that H'Q) > 0, we conclude that H is a nonnegative 
function. Let r\ v be given by the sequence [12] 

r \x\< v- % 

%{x) = I log (v 2 \x\) I log (v 2 R) v~ 2 < \x\ < R 
{ 1 R<\x\. 

Notice that r\ v vanishes on a neighborhood of the origin for any finite parameter 
v; but as v tends to infinity, r\ v converges pointwise to 1, whereas Vfj v converges 
to zero in L n {D). 

Inequality (189|) now assumes the form 

0< / [(£ + f)H]-( V r]„) 2 H*l 
Jd 

= f {[d k (l + a3 k )d J+ f]H}{ mv ) 2 H*l = 
Jd 

[1 + a ]k ) djH ■ d k \{m v f H] * 1 + / fH 2 (rfrj v f * 1, 

(i + o»**) djH ■ d k \{rfnuf h}*\ = 



where 



D 
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-2 / (l + a> k )Tjlr](dkV)H -djH*\ 

J D 

-2 f (l + a^^rjAdkV^H-djH*! 
Jd 

- [ (1 + a jk ) (r,r} u ) 2 {d 3 H ■ d k H) * 1 
= - (2i x + 2*2 + i 3 ) . 

That is, 

k = [ (1 + a jfe ) (wJ 2 ■ * 1 

<2(|i 1 | + |i a |)+ / (m v ffH 2 *i. 

Estimating the integrals on the right individually, we have 

2*i = 2 / (1 + a J ' fc ) r^jj (^) H-djH*l< 
Jd 

e f (v^f \WH\ 2 * 1+ 
Jd 

- I f v \Vr,\ 2 H 2 *l. 
£ Jd 

2^ 2 - 2 f (1 + ^' fe ) ry 2 ^ (d k rj v ) H-djH*l< 
Jd 

el M) a |Vi/| 2 * 1 + - f V 2 \V7]„\ 2 H 2 *1 
Jd £ Jd 

1 

= £«21 + -«22, 

£ 

where the small constants £ depend on the constant C in the upper inequality 
of (HE]), and 

* 22 =/ ?? 2 |Vr7j 2 i/ 2 *l<C||V77j| 2 ||r,i7|| 2 n/(ll _ 2) . (90) 

The right-hand side of (f9"0"| tends to zero as ^ tends to infinity. Absorbing small 
constants on the left, now using the lower inequality of (55]) , we conclude that 

(1 - 2e) f (r^) 2 |Vif| 2 * 1 < 
Jd 

C{e) Urj 2 u \Vr)\ 2 H 2 * 1 + \\Vrj v \\l ll^||L/(„-2) + ll/IU/2llw^llL/ (n -2)) » (91) 
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where 



C 



\\VV,H\\in/(n-2) < C J d |V (vVuH) | 2 * 1 < 

f / IV (W.) \ 2 H 2 * 1 + / (wJ 2 |Vif| 2 * l) 



(92) 




The first integral on the extreme right-hand side of inequality ([92]) has essentially 
already been estimated in (|90|) . and the second can be subtracted from the left- 
hand side of (f9"Tj) provided that its coefficient in (f9"Tj) . the L"/ 2 -norm of / over 
f2, is sufficiently small. Letting v tend to infinity, we have in the limit 



We conclude that H is a weak solution in a neighborhood of the singularity. 
This completes the proof of Theorem 6. 

7.3 Proof of Theorem 8 

1. Outline: The result follows from the application of Lemmas 3 and 4 and 
Theorem 5, above, to the proofs of Theorem 6 and Corollary 7 of [27] . The 
absence of condition (|3"3"|) results in a small change of L p conditions on the 



2. Details: Initially proceed as Sec. 7.2, but choose the test of functions r\ 
to be a sequence of functions possessing the following properties: 

(a) 7y„ e [0,1] Vi/; 

(b) ^ = in a neighborhood of E W; 

(c) lim^oo T) v = 1 a.e., 

(d) lim™||V^J| £ n- m - e =0 

(c.f. [33], Lemma 2 and p. 73). The function T is given by [32] 




solution. 



T{H) = 



(l/7i) [72^ 2 " 71 ^ 71 - (7i " 72) ^ 72 ] 



< H < t, 
H>1. 



The functions T{H) and 



G(J?) = T(H)F'{H) - 72 



satisfy ([32], p. 



280; see also Sec. 3 of [12]): 



31 



\G(H)\ < T(H)T'(H); 
G'(H) > CT'iHf. 
Replace the test function in Sec. 7.2 by the test function 

where rj is defined as in Sec. 7.2, and substitute this value into ([55)) . We obtain 



or 



0< / {[d k (l + a? k )d j + f]H}{r^ v fG{H)*l = 

JD 

f (1 + a? k ) djH ■ d k \( V f} u ) 2 G(H)} * 1 + / fifrfvH ■ G(H) * 1, 

JD 1 J JD 

(1 + a? k ) d k (rflj v ) (djH) G(H) * 1+ 
/ (l + a^irfrivfG'iH^HdkH + lK [ f(r,r} v ) 2 H-G(H)*l. 

JD JD 



I D 

Writing this inequality in the short-hand form 



h+h<h, (93) 

we proceed analogously to inequalities (28)-(35) of [26] , making the following 
changes from the notation of [55] to our notation: ifj — > rj, Q — > H, H — > T , 
S — > G, $ — * /. Explicitly, 

h = 2 / (1 + a> k ) {rfij v ) 8 k {m v ) ■ G{H)8 H * 1 

JD 

>-C I {rjfjv) IV (mu)\\HH)VHH)\ * 1 

JD 

> s [ im,) 2 \^HH)\ 2 * 1 - C(e) [ |V {rflf) \ 2 T{H) 2 * 1; 

JD J D 

h>C f {m„?\F'{H)\ 2 \VH\ 2 *l 

JD 

>C [ (wj 2 |VJ^(ff)| 2 * 1; 

JD 

h = I f(r,r)u) 2 H-G(H)*l< 

JD 

{ \f\{m v f\H-T'{H)\\T{H)\*l< 

JD 

72 J d |/| (r/r/J 2 \HH)\ 2 * 1 < ||/|| n/2 Qf |w^(ff)| 



2n/(n-2) + x 



(n-2)/n 

(94) 
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As in ([9"2"|) , we apply the Sobolev inequality to the right-hand side of this ex- 
pression, followed by the Minkowski and Schwartz inequalities: 



< C 



( / \m„HH)\ 2nKn - 2) * ij < cj D |v [w^(if)] | 2 * i 

/ \(W V )r]^(H)\ 2 *l+ f \(Vrj v )r)F(H)\ 2 *l+ [ {m„?\VF{H)\ 2 *l 

JD JD JD 



hi + I32 + h 



The term ^33 can be subtracted from the left-hand-side of inequality ([93]) , as its 
coefficient in (IM1) . the L"/ 2 -norm of /, is small on small discs as a consequence 
of condition (1561). 



*32 



|(V^)^(i7)| 2 *l<C( 7 i,72^) / KV^f^H 2 ^ 

JD 



* 1 



<C\\Vrj v \\ 2 n _ m _ e \\H 2ll \b- 

Letting v tend to infinity, the term on the right-hand side is zero for every value 
of I. Now letting I tend to infinity and using Fatou's inequality, we conclude 
that 

/ tfl^iH^l 2 *1<C f \Vr)\ 2 H 2 ^*l. 

JD JD 

Apply Theorem A.l, taking U = H~< 2 . Then U S H l > 2 (D) and W — U x satisfies 
inequality ([87]) for A = 1/72- Because 72 > 1/2, we can conclude that A < 2 as 
required by Theorem A.l. We now want to check that we can choose 72 < 1, 
in order to obtain A > 1 as also required by Theorem A.l. Because H lies in 
the space L 2071 (L>) n L 2 ^ 2 (D), we let 7i/3 = 72. Substituting the definition of 
[3, we find that we can choose 72 < 1 for 71 > 1/2 provided m + 4 < n, which 
is satisfied by hypothesis. This completes the proof of Theorem 8. 
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